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The distance modulus determined from Carmeli's cosmology fits 
the accelerating universe data of the high-redshift type la 
supernovae without dark matter 

John G. HartnetlQ 

School of Physics, The University of Western Australia 
35 Stirling Hwy, Crawley 6009 WA Australia 
(Dated: February 2, 2008) 

The velocity of the Hubble expansion has been added to General Relativity by Moshe Carmeli 
and this resulted in new equations of motion for the expanding universe. For the first time the 
observational magnitude-redshift data derived from the high-z supernova teams has been analysed 
in the framework of the Carmeli theory and the fit to that theory is achieved without the inclusion 
of any dark matter. Best fits to the data yield an averaged matter density for the universe at 
the present epoch Q, m ~ 0.021, which falls well within the measured values of the baryonic matter 
density. And the best estimate of Qa + O m « 1.021 at the present epoch. The analysis also clearly 
distinguishes that the Hubble expansion of the universe is speed-limited. 

PACS numbers: 95.30.Sf 95.35.+d 98.62.Py 98.80.Es 



I. INTRODUCTION 

The metric 0, Q used by Carmeli in a generally covari- 
ant theory (Cosmological General Relativity) extends the 
number of dimensions of the universe by the addition of a 
new dimension - the radial velocity of the galaxies in the 
Hubble flow. The Hubble law is assumed as a fundamen- 
tal axiom for the universe and the galaxies are distributed 
accordingly. 

In determining the large scale structure of the universe 
the usual time dimension is neglected as observations are 
taken over such a short time period compared to the mo- 
tion of the galaxies in the expansion. In this case, the 
usual time coordinate (t) does not need to be considered. 
This leaves only four dimensions to be considered - three 
of space and one of velocity. In general, however, the 
new 5D cosmology developed by Carmeli contains all of 
general relativity as a subset. All of the results in gen- 
eral relativity that are experimentally supported are also 
obtained in Cosmological General Relativity (CGR) 

When discussing the motions of stars and gases in 
galaxies the five dimensional cosmology must be used. 
Both the usual geodesic equation and a new phase 
space equation have been derived using CGR. From this 
Carmeli was able to derive a Tully-Fisher like equation 
3] , which indicates that the existence of halo dark matter 
is not necessary to be assumed in spiral galaxies. 

It is then worth seeing to what extent the Carmeli 
cosmological theory fits the observed data. This paper 
fits the theory to the observed magnitude-redshift data 
taken from the high-redshift type la supernovae observa- 
tion teams. From this analysis on the largest scales of the 
universe no dark matter is necessary when the Carmeli 
model is assumed. 



A. Cosmological General Relativity 

Here the CGR theory is considered using a Riemannian 
four-dimensional presentation of gravitation in which the 
coordinates are those of Hubble, i.e. distance and veloc- 
ity, or more precisely, proper distances as measured by 
the Hubble law and the measured redshifts of galaxies. 
This results in a phase space equation where the observ- 
ables are redshift and distance. The latter may be de- 
termined from the high-redshift type la supernovae (SNe 
la) observations. 



B. Line element 



In the case considered here the line element is 



ds 2 = r 2 dv 2 



e ( dr 2 -R 2 (d0 2 + sin 2 Bdcb 2 ) ( 1 ) 



'Electronic address: john@physics.uwa.edu.au 



where dr 2 = (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 , £ and R are 
functions of v and r alone and comoving co-ordinates 
x^ = (x , x , x 2 , x 3 ) — (Tv,r,6,4>) are used. The new 
dimension (v) is the radial velocity of the galaxies in 
the expanding universe, and is not the time derivative of 
the distance. The parameter r, the Hubble-Carmeli time 
constant, is a constant for all observers at the same epoch 
and its reciprocal (designated h) is the Hubble 'constant' 
measured in the limit of zero gravity and zero redshift, 
which is only approximately the Hubble constant Hq. 

Equation then represents a curved spacevelocity 
which, like in general relativity, may be represented by 
a four-dimensional Riemannian manifold with a metric 
g^u and a line element ds 2 — g^dx^dx" . This differs 
from general relativity in that here the x° coordinate is 
velocity-like instead of time-like as is the case of x° = ct, 
where c is the speed of light and a universal constant. 
The parameter t is the time coordinate. In this theory 
x° = tv, where r is also a universal constant. The other 
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1,2,3, are space-like, as in gen- 



three coordinates x k , k 
eral relativity. 

The line element represents a spherically symmetric 
isotropic but not necessarily homogeneous universe. The 
observations are made at a definite time and therefore 
dt = does not appear in QJ. 



C. Field equations 

In CGR, as in general relativity, one equates geometry 
to physics. In this theory Einstein's field equations 



(2) 



are modified. 

The energy-momentum tensor (T^) takes on a differ- 
ent physical meaning. The coupling constant (re) that 
relates the geometry G M „ in J5J to the energy terms T^ v 
is also different. However the form of (J2J is exactly the 
same as in general relativity, but with re = Snk/r and 
k = Gt 2 /c 2 where G is Newton's gravitational constant. 
Therefore re = 8itG/c 2 t 2 . 

The correspondence with general relativity is easily 
seen by the substitutions c — * t and t — > v. In this new 
theory the energy-momentum tensor T^ v is constructed 
with these substitutions. As usual T^ v — pu^u 1 " ', where 
where p is the average mass/energy density of the uni- 
verse and m m = dx^ /ds is the four-velocity. 

In general relativity Tq = p. In Newtonian gravity 
the potential function is defined by the Poisson equation 
V 2 = AnGp. Where p = the vacuum Einstein field 
equations are usually solved in general relativity and sim- 
ilarly Laplace's equation in the Newtonian theory. 

These are valid solutions but in cosmology there never 
exists the situation where the density p is zero because 
the universe always contains matter and energy. So in 
order to equate the rhs of J3J) to zero Carmeli took Tq ^ 
p but Tq = p e ff = p — p c , in the appropriate units. 
Here p c is the critical or "closure" density and in this 
model p c — 3/8itGt 2 w 10~ 29 g.cm~ 3 . Therefore in 
CGRF"=p e// iiV. 

The result is that we can view the universe, in spaceve- 
locity, or phase space as being stress free when the mat- 
ter density of the universe is equal to the critical den- 
sity. That is, the effective density p e ff — 0. This then 
gives us the analogous situation to that in the Newto- 
nian and Einsteinian theories. Besides the assumption 
of the universality of the Hubble Law, this is the second 
fundamental assumption in this cosmology. 



D. Phase space equation 

In spacevelocity the null condition ds = describes 
the expansion of the universe. In the limit of no mat- 
ter or gravity (i.e. = 1) the null condition yields 



dr/dv = ±t, which when integrated with appropriate ini- 
tial conditions results in r — tv in an expanding universe. 
This can be rewritten as v = hr rs H$r, the Hubble law 
in the zero-gravity limit. 

It follows from Q for a spherically symmetric isotropic 
distribution of matter where spherical spatial coordinates 
(r, 9, </>) are used and taking into account d6 — d<j> = 
(the isotropy condition) that 



r 2 dv 2 - e s dr 2 = 



which results in 



dr 
dv 



±re'^ 2 . 



(3) 



(4) 



The positive sign is chosen for an expanding universe, 
which was solved in 0, [j| . 



E. Solution to field equations 

Carmeli found a solution to the resulting field equa- 
tions 3, with the necessary condition that R' > 0, as 



and 



R = r 



1 



(5) 



1 + f(r) 



(6) 



where f(r) is an arbirary function of r and satisfies f(r) 
1 > 0. The solution is where 



f(r) 



i - n 



(7) 



and Q = p/p c . The matter density f2 is assumed to be the 
smoothed average density for matter that is evenly dis- 
tributed throughout the universe. However because we 
look back through past epochs the density is a function 
of the redshift z. 

By substituting © with Q into Q we get Carmeli's 
result 



dv 

_=ri/l + (l-n). 



(8) 



where the positive solution has been chosen for an ex- 
panding universe. 

Equation ||SJ| may be integrated exactly to get 



r(v) 



sinh -VT^n Vfi. (9) 



Thus © may be written in terms of normalized or nat- 
ural units r/cr and for arbitrary z = v/c, 



r sin 



(10) 
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Considering the expansion of the universe it is clear that 
HI U| l describes a tri-phase expansion. Initially the uni- 
verse is very dense and 51 > 1 so the hyperbolic sine func- 
tion becomes a normal trignometric sine function describ- 
ing a decelerating expansion. Then the density ft reaches 
unity and the rhs of (|10|l becomes equal to z which de- 
scribes a coasting stage. Finally the density decreases 
and fl < 1 as the universe continues to expand. Then 
the hyperbolic sine function represents an exponentially 
accelerating universe. 

There are four symbols used in this paper for density 
expressed as a fraction of the critical density. The sym- 
bol fl represents the matter density at any epoch defined 
by redshift z and because we take dt = 0, fl is there- 
fore only a function of z. The symbol fl m represents 
the matter density at the present epoch and therefore is 
a constant on the time scale of any measurements used 
here. The symbol fib specifically represents the baryonic 
matter density at the present epoch, which in this pa- 
per we show is identical with fl m . The symbol fl\ is the 
vacuum or 'dark' energy contribution to gravity and is a 
function of redshift z. In fact, because the cosmological 
constant (A) does not appear explicitly in Carmeli's cos- 
mology this latter component is really a property of the 
metric. 

The structure of this paper is as follows: SectionlTlldis- 
cusses the variation of matter density with redshift and 
presents an equation, modified by the author [TT| from 
the Carmeli equation, that describes proper distance as a 
function of redshift in the Carmeli cosmology. Section lTTll 
is where the new work begins. In this section, for the first 
time, a magnitude-redshift relation is fitted to the data 
taken from the high-z supernovae teams for 
a model that assumes flat spacevelocity. In Section IIVI a 
magnitude-redshift relation is fitted to the same data but 
this time for a model that assumes curved spacevelocity. 
SectionlVlcalculates the value of h = 1/r. Section lvTl cal- 
culates the value of fl^ and the total £1^ + 51 as a function 
of redshift. The appendix calculates, in a more rigorous 
way than in Section llVl the effect of curved spacevelocity 
on matter density as a function of redshift. 



II. MATTER DENSITY VERSES REDSHIFT 

Now let us consider the density of matter as a function 
of redshift, z. Carmeli assumed that the value of fl in 
HI U| l is fixed and plotted curves as functions of redshift 
for various values of fl. See figure A4, on page 134 of 0. 

In 1996 he predicted that the universe must be acceler- 
ating and simulated the form of the high redshift data 
of Riess et al 01, published in 1998, which announced an 
accelerating universe following the observations of Gar- 
navich et al [t| and Perlmutter et al [l5| . 

Therefore in 1998, Carmeli assumed a value of total 
matter (normal bayonic matter + dark matter) density 
fl = 0.245, which was the accepted value then for fl m 
the matter density at the present epoch. However, more 



correctly fl varies as a function of redshift, z. 

Carmeli never fitted his theory to any of the high red- 
shift SNe data with a least squares method or any other, 
nor did he determine a value for the matter density fl m 
himself. 



A. Baryonic matter density 

The density of normal baryonic matter in the universe 
has generally come from considerations of big bang nu- 
cleosynthesis (BBN)^I, [^|. The predicted primordial 
abundance of the light elements are assumed to have been 
produced in the big bang and not in stars. Hence a com- 
parison between predicted abundances and abundances 
inferred from observed baryon-to-photon ratios are used 
to put an upper limit on the baryon density. 

Deuterium is believed to be only destroyed by stellar 
processing, therefore any observation of its interstellar or 
solar system abundance would put a lower bound on its 
primordial abundance since it is believed that it was all 
created in the BBN process. As a result a lower bound 
on deuterium translates into an upper bound on the den- 
sity of baryons. Following this line of logic the observed 
interstellar value of D/H > (1.6 ± 0.2) x 10-^_puts a 
limit on the baryon density fl b h 2 < 0.027(2<r) 0- An 
earlier determination produced flbh 2 < 0.015±0.005 [l4|. 
(Here fib is the baryonic density expressed as a fraction 
of the critical density and h is the Hubble constant as 
a fraction of 100 fcm.s _1 Afpc _1 and not to be confused 
with h = 1/r used in this paper.) Assuming a value 
of the Hubble parameter h = 0.7 these two sources put 
fib < 0.055 and < 0.031 respectively. 

The aforementioned values of fib agree well with the 
locally measured baryonic budget. One study yields a 
range 0.007 < fib < 0.041 at z w with a best guess of 
fib ~ 0.021 where a Hubble constant of 70 km.s -1 Mpc -1 
was assumed @- 

Now, determinations of the total matter density (in- 
cluding dark matter), at high redshifts, have been made, 
for example, from X-ray measurements of the gas content 
in g alaxy clusters combined with the total virial masses 
[2Cj . The calculation of the virial masses in turn assumes 
that the equations describing the motions of the con- 
stituents in the clusters are correct and therefore gravita- 
tional potentials are assumed that generate masses larger 
than the baryonic matter content. Hence a significant 
dark matter content results. 

However it is well-known that mass-to-light ratios vary 
over all mass scales from galaxies to superclusters [21|. 
This may be the result of an incorrect understanding 
of the physics which has resulted in incorrect dynami- 
cal masses. That conclusion is suggested by Carmeli's 
derivation of a Tully-Fisher type relation for galaxies Q 
and by a new post-Newtonian equation describing the 
dynamics of the stars and gases in spiral galaxies |10| . 
The latter produces the galaxy rotation curves, without 
any dark matter, where previously it was necessary to 
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assume that they are evidence for halo dark matter. 

Therefore the application of the Carmeli theory to the 
very important problem of galaxy cluster dynamics is 
necessary to see if it also eliminates dark matter on that 
scale. Eventually it is necessary that the Carmeli the- 
ory also determine cluster masses from the virial theo- 
rem, observed X-ray temperatures, and from weak gravi- 
tational lensing with results consistent with each method 
and with the observed baryonic matter density. However 
such an exercise is beyond the scope of this paper. 

The question that remains for this paper to answer is 
"Do the high-z SNe la data ht this model without the 
need to assume the existence of dark matter?" In future 
research, other issues can be addressed, for example, to 
explain the much higher than expected dynamical masses 
of clusters at high redshift, and the shape of the black- 
body spectrum of the cosmic microwave background as 
well as the acoustic peaks in its spatial spectrum. 

B. Matter density at the present epoch 

This paper makes no determinations about the mat- 
ter density, dark, baryonic or otherwise at high redshift 
except the following, with a justification in section IVTI 

n = a m (i + z) 3 , (ii) 

where f2 m is then the averaged matter density, at the 
present epoch, expressed as a fraction of the critical or 
"closure" density. f2 m can be considered to be constant 
on the time scale of any measurements used here. Equa- 
tion l|ll|) results from the fact that as the redshift in- 
creases the volume changes as (1 -I- z)~ 3 assuming flat 
Euclidean space. 

In the standard Friedmann-Lemaitre (F-L) theory the 
equivalent expression is 

Q = n m a- 3 /H(a) 2 , (12) 

where a = (1 + z)^ 1 the scale factor (because at the 
present epoch, as usual, a is assumed to be unity) and 
H{a) is the Hubble term which quantifies the curvature 
of the expansion and defined H 2 = (a/a) 2 , where the dot 
is the time derivative. 

In this theory the curvature results from the inclu- 
sion of the new velocity dimension and the implicit as- 
sumption of the universality of the Hubble law. This 
introduces curvature through the concept of spaceveloc- 
ity, which is discussed later. (See section IIVI and the 
appendix. The time derivative of a scale factor is not 
relevant.) 

In the first instance I assume a flat spacevelocity model 
in (jnj, which is identical to JTJJ where H(a) = 1. This 
would correspond to a dust dominated spatially flat uni- 
verse in the standard theory. 

Now the value of the baryonic matter density com- 
monly cited is p — 3 x 10~ 31 g.cmT 3 . If this is assumed 
to be the total matter density then fi m w 0.03 in JUJ. 



Substituting ((TIT) mto G3 we S et 

sinh (zy/l - n m (l + z) 3 ) 

— = \ =—*-. (13) 

v /l-O m (l + z)3 

It has been shown pl | that dm with n = 0.245 (which 
Carmeli initially assumed) and 113fl with fi m = 0.03 are 
practically identical over the redshift range < z < 1. 
The difference between the two equations over the do- 
main of the measurements is much less significant than 
the fit to the data. 

Once it is shown that this theory fits the actual ob- 
servational data it follows that this effectively eliminates 
the need for the existence of dark matter on the cosmic 
scale. 



III. COMPARISON WITH HIGH-Z TYPE IA 
SUPERNOVAE DATA 

In order to compare (|13|) with the data from the high 
redshift SNe la teams the proper distance is converted to 
magnitude as follows. 

m(z) =M + 5logV L {z;fl), (14) 

where T>l is the "Hubble constant free" luminosity dis- 
tance, is related by 

V L (z;Sl) = (l+z) — (15) 

CT 

using CS|| and 

M = 5log(cr) + 25 + M B + a. (16) 

Here M. incorporates the various parameters that are 
independent of the redshift, z. The parameter Mb is 
the absolute magnitude of the supernova at the peak of 
its light-curve and the parameter a allows for any un- 
compensated extinction or offset in the mean of abso- 
lute magnitudes or an arbitrary zero point. The absolute 
magnitude then acts as a "standard candle" from which 
the luminosity and hence distance can be estimated. 

The value of Mb need not be known, neither any other 
component in M., as M. has the effect of merely shifting 
the fit curve (|15(l along the vertical magnitude axis. 

By choosing the value of r = 4.2 x 10 17 s — 
13.3 Gyr, which is the reciprocal of the chosen value 
of the Hubble constant in the gravity free limit h = 
73.54 km.s' 1 Mpc~ 1 (see sectionE} M = 43.06+Ms+a. 
In practice, where the distance modulus (to— Ms) is used 
in the curve fits, a is an arbitrary free parameter. 

It has been one of the goals of cosmology to determine 
if the Hubble expansion is speed limited. That is, to an- 
swer the question of whether the cosmological expansion 
is governed by the relativistic Doppler effect or not. If it 
isn't then in 1|13|) v/c = z (as written), but if it is then 
as v approaches c it is necessary to replace z in l(T3|) with 
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FIG. 1: (Color online) Data taken from Table 5 of Riess et 
al. Curve 1 (purple solid line) is the best fit and respresents 
the speed limited model. Curve 2 (orange dots) is the non 
speed limited model. Both curves use the same value of Q m — 
0.02. In order to make curve 2 pass through the highest data 
points the required mass density Q. m must be larger than 0.10. 
However such values are unphysical 



v/c = ((1 + z) 2 - 1)/((1 + z) 2 + 1). So by fitting to the 
data of the high-z SNe la it is possible to test this ex- 
perimentally. Behar and Carmeli 1] have asserted that 
latter (speed-limited version) is the correct form, so this 
then becomes an experimental test of the prediction. 

In figure ^ GH has been plotted against the distance 
modulus (m — Mb) data taken from Table 5 of Riess et 
al [l||. Curve 1 is the best fit to the data using 
with v/c = ((l + z) 2 -l)/((l + z) 2 + l), the speed-limited 
model and curve 2 is 113fl with v/c = z the non-speed- 
limited model. This immediately answers the question. 
The Hubble expansion, as seen from our position as ob- 
server, is speed-limited, by the usual relativistic Doppler 
effect, which confirms the theoretical prediction. 

It can be seen that curve 2 departs from curve 1 slightly 
above redshift, z — 0.5. Both curves have been plotted 
with a value of Q, m = 0.02, which is the best fit for curve 
1. However different 'best fit' offsets (a) have resulted. 
For curve 1, a = 0.41 and for curve 2, a = 0.011. There 
are a total of N — 185 data and the residual for the speed- 
limited version of U3J) is R = 0.9960, with x 2 = 19.94 and 
hence x 2 /N — 0.1078. The residual for the non-speed- 
limited version of U3 is R = 0.9927, with X 2 = 37.465 
and hence x 2 /N = 0.2036. 

In order to make curve 2 pass through the highest data 
points the requirement on the mass density is fl m > 0.10, 
which is unphysical assuming only baryonic matter. As 
mentioned in section^ the local baryonic matter budget 
has been measured in the range 0.007 < Q& < 0.041 Q 
with a best guess of fib ~ 0.021. The latter is consistent 
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FIG. 2: (Color online) Data taken from Riess et al for z > 0.2. 
Curve 1 (black solid line) is the best fit and respresents I14H 
with the speed-limited 117H . Curve 2 (orange dots) represents 
114H with Q = fi m = and curve 3 (purple dot dash) repre- 
sents 1141 with Q, — 1 



with the best fit from this data set. Therefore, in the 
following analysis, £l m = Qb i-e. no dark matter needs 
to be included. Also only the speed-limited form of 113|) 
will be used. That is 



sinh 



^i-n m (i + z)^ 

^/l-{l m (l + zf 



(17) 



where ? = ((1 + z) 2 - 1)/((1 + zf + 1). 

Figure |2 shows the data of figure ^ but for z > 0.2. 
Curve 1 (solid line) is the best fit curve over the range of 
the selected data and also consistent with measured bary- 
onic density i.e. fi m > 0.007. Instead of fi m = 0.02, the 
best fit curve requires fi m = 0.007 ± 0.050 (statistical). 
However, the value of f2 m may range up to 0.04 with little 
change in residuals. A slightly increased value of a = 0.53 
results compared with the fit in figure ^ This may indi- 
cate a small 0.12 additional extinction over the low red- 
shift data. Here N = 106, the residual R = 0.9946, 
X 2 = 17.107 and hence X 2 /N = 0.1614. Curve 2 (dots) 
represents ifEI}) with £1 = Q m = and curve 3 (dot dash) 
represents (|14[1 with fi = 1. Curve 3 represents a universe 
where the matter density is always critical regardless of 
epoch. 

Figure [3] shows the 8 new SNe la data from Tonry et 
al 19] for sources with z > 0.2. The resulting best fit 
is shown, with fi m = 0.073 ± 0.098 (statistical), a — 
0.56, the residual R = 0.9938, x 2 = 0.1460 and hence 
X 2 /N = 0.0183. Tonry et al used 4 methods of analysis 
and averaged the results. The data shown in figure [3] is 
drawn from their Table 6, column 4. 

Figure0|shows 30 type SNe la data (which includes 11 



6 



44 



CO 
T3 

o 



CD 42 
O 

c 

Cfl 



co 



40 




0.2 



0.4 



0.6 0.8 

Redshift (z) 



1.2 



FIG. 3: Data taken from Tonry et al. The curve (solid line) 
is the best fit and respresents (f 1 4^ with the speed-limited ff 1 7|l 
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FIG. 4: Data taken from Knop et al for z > 0.2. The curve 
(solid line) is the best fit and respresents 11411 with the speed- 
limited JHJ 



new SNe) from Knop et al ^3 f° r sources with z > 0.2. 
In this case the magnitude has been plotted against red- 
shift. These data are taken from column 4 of Tables 3 
and 4 in Knop et al. As with the other data various 
corrections have been applied, including i\-correction, 
Galactic extinction, and light-curve stretch correction. 
(See references for details). In figure 0] the host galaxy 
corrected data (column 5 from (Knop et al)) was not 
used as it added a large scatter without significant im- 
provement to the fit. The resulting best fit is shown, with 
n m = 0.019±0.137 (statistical), the residual R = 0.9217, 



X 2 = 1.9323 and hence x 2 / N = 0.0644. Because the ef- 
fective apparent magnitude mseff was used instead of 
distance modulus the offset is Mb + a = —18.760. 

By taking the statistically weighted average of the mat- 
ter density fi m obtained from each curve fit in figures |3 
- 01 we get an average value 



a, 



E 



7 



i 



0.021 ±0.042, 



(18) 



where fi m i are the best fit values from the three data sets 
used and a is the standard error for fi m i. 



IV. CURVED SPACEVELOCITY 

In section|n]on matter density it was assumed that the 
density could be described by the function in ijTTJl that 
assumes the universe is essentially spatially Euclidean. 
This assumption is justified in section IVII However it 
is obvious that besides for small z and where fi = 1, 
spacevelocity is not flat in general and is expected to mod- 
ify the expression in jH) . When fi > 1 it is curved and 
closed and when 17 < 1 it is curved and open. Equa- 
tion (jTT|) relates the matter density at any epoch to 
the present epoch value (Q m ) in Euclidean space or flat 
spacevelocity. 

Rewriting in an expanding universe we get 



cr dz 



(19) 



The result is the gradient of the spatial co-ordinate r 
with respect to the redshift z, normalized by the Hubble 
length ct. 

Therefore it is clear from l|19|) that e^ 2 = 1 under the 
conditions that produce the Hubble law in the limit of 
weak gravity or for flat spacevelocity. It can be seen from 
© and that e? /2 = 1 when = 1. Therefore we 
can use e^ 2 to define the curvature of spacevelocity as 
a function of redshift z. Because the density scales as 
the inverse cube of the radial coordinate and since the 
radial coordinate scales differentially as e - ^/ 2 , then e 3 ^/ 2 
may be a good estimate of how the density scales with 
redshift z. 

See the appendix for a more rigorous approach. The 
more rigorous approach yields an integral that can only 
be approximated and results in a transcendental equation 
which cannot be used in the curve fits. Whereas the 
former approach results in an analytic function that can 
be used in the curve fits. 

By combining i|ll[l with |jfj| and J7J|, for e^' 2 ^ 1 we 
can write 



n m (i + z) 3 



(l + (l_0)(£) 



3/2 ' 



(20) 



When = 1 in the denominator of l|20|l we recover l|ll|l . 
To get the density fi as a function of z we substitute r j ct 



7 



0.8 



G 



% 0.6 
c 

CD 
T3 

^ 0.4 

"5 



0.2 

















q 
o 

/l 
















,/ / 







.- 












i Lb i ■ 


. - 




R 



0.5 1 1.5 

Redshift, z 



-18.76 and residual i? = 0.9218, \ 2 = 1-9317 and hence 
X 2 /N = 0.0644. 

Finally by taking the weighted average as before 

fim = V— ^/V -4 = 0.020 ± 0.068. (22) 
erf ^ erf 

The new density function 1(21(1 yields the same result 
as the previous model. However in the Appendix it is 
shown that a more rigorous approach indicates that the 
curved spacevelocity model matter density approximates 
the density of the flat spacevelocity model with only a 
1.28 multiplying factor for z < 2. This means the best 
fit derived value Vt m = 1.28 x 0.021 = 0.027. However 
the difference in the fits using Q, m = 0.021 or 0.027 is 
extremely small and only apparent for the few data in 
the region z > 1, but much less than their associated 
error bars. 



FIG. 5: (Color online) Matter density as a function of redshift 
z for the two models. Curve 1 (black dots) is the flat spacevel- 
ocity density from lilt with Q m — 0.021 taken from list . The 
curved spacevelocity density from 12 U with Q m = 0.021 is rep- 
resented by curve 2 (orange dot dash), and with f2 m = 0.027 
by curve 3 (purple solid line) 



from 1(17)1: 

n = 



n m (i + zf 



(l + smh? (sy/l -f! m (l + z) 3 )) 



3/2 ' 



(21) 



where <;=((! + zf - 1)/((1 + zf + 1). 

This results in a new density function which is com- 
pared with (|11|) in figure^ Curve 1 is the flat spaceveloc- 
ity density from 1(11(1 with fi m = 0.021 taken from (|18fl . 
Curve 2 is for the curved spacevelocity density from 1(21(1 
with the same value of fi m = 0.021. Note the curves 1 
and 2 are equivalent for z < 0.4. For small z we recover 
l|H)l from J2D. Because of the condition f(r) + 1 > 
equation 121|) is valid for all z < 4 provided Q, m < 0.03. 

The density from l|21|) is then substituted into i|17|) 
replacing ft m (l + z) to get a new expresssion for r/cr 
with only one free parameter Q m . This was then used in 
the same curve fits in section IIHI (figures ^ - EJ with the 
following results. 

In each case the curve fits were improved by the new 
density model except for the figure |21 fit, which had a 
marginally worse y 2 statistic. For the curve 1 of figure 
n the fl m = 0.02 fit had a reduced but the same R = 
0.9960. This resulted in a slight change in the offset a = 
0.41. For figureHthe best fit yielded O m = 0.016 ±0.115 
(statistical) with offset a = 0.56 and residual R = 0.9945, 
X 2 = 17.208 and hence X 2 /N = 0.1623. For figureElthe 
best fit yielded TO = 0.020 ± 0.097 (statistical) with 
offset a = 0.56 and residual R = 0.9928, \ 2 = 0.1684 and 
hence x 2 / N = 0.0211. For figure 01 the best fit yielded 
fl m = 0.029 ± 0.175 (statistical) with offset M B + a = 



V. HUBBLE PARAMETER 

By inverting i|17[l , multiplying both sides of the result- 
ing equation by v j c and using the expression v = Hqt we 
get the following 



Hn = h 



sinh (<;VT~~n) 



(23) 



with the density taken from either (jlljl or 1)21(1 . 

As a result we have an expression which indicates that 
the Hubble parameter H is scale length related and a 
function of the universal Hubble constant (h). It has 
been common to find figures cited in the literature that 
indicate 'long' and 'short' scales for Ho. Studies with the 
Hubble Space Telescope using classical Cepheid variables 
yielded H = 80 ± 17 km.s~ 1 Mpc~ 1 whereas studies 
usingSN la yielded H = 67 ± 7 kms^Mpc' 1 0. See 
also @. 

It then follows from that h w 73.54 km.s^ 1 Mpc~ 1 
for fi m = 0.021 and the flat spacevelocity model i|ll|) 
when a value of Hq = 70.00 km.s^ 1 Mpc^ 1 at z = 1 is 
assumed. In any case, the specific value of h or it re- 
ciprocal t cannot be determined fom the data fit. An 
independent method must be used to determine the pre- 
cise value of t. 



VI. DARK ENERGY 

The vacuum or so-called 'dark' energy parameter S1a 
does not appear explicitly in Carmeli's CGR. Hence the 
term 'dark' energy is a misnomer, and probably 'vacuum' 
energy is more correct. In any case, it is really a property 
of the metric. Only by a comparison with the standard 
F-L models can an assignment be made P, Q . 

The vacuum energy density pa = A/8irG (in CGR) = 
3Hq/8itG (from the standard theory). Also in CGR the 
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FIG. 6: (Color online) Qa (curve 1), Q (curve 2) and the total 
mass/energy density f^A + fi (curve 3) as a function of redshift 
2 for the flat spacevelocity density model with fl m — 0.021 

critical density p c — 3h 2 /8irG. Therefore fl\ = pa/ Pc = 
(H /h) 2 and it follows from that 

where the density fl is taken from and the resulting 
Qa, as a function of z, is shown in figure for the flat 
spacevelocity density model. Curve 1 and 2 are respec- 
tively the values of Oa and f2 with f2 m = 0.021. 

Curve 3 in fig. shows the values for the total energy 
density Q + Q\ as a function of redshift, z, for the flat 
spacevelocity model with Q m — 0.021. The total density 
fl + Q\ ps 1 for z < 1 and still remains close to unity 
up to z = 2. This means out to z w 2 the universe is 
quasi-Euclidean. 

Table [I] shows the two models compared with D, m = 
0.02,0.03,0.04 which are all within the bounds of the 
measured baryonic matter density at the present epoch 
(z w 0). Fits from the two models (flat with fl m w 0.02; 
curved with fi m w 0.03) both result in approximately the 
same total density f2 + Ha = 1.067 at z = 1. 

For small z the total density becomes 

a + n A » (i + fi m ) + 3^o m . (25) 

It follows from that for Q m w 0.021 at z w the 
total density fl + Q A « 1.021. 

Furthermore now let us consider the time development 
of these densities in a qualitative sense only. From 124|) 
it follows that as the universe expands the total density 
tends to the vacuum energy density J7a ~~ * 1 because 
— ► 0. This means a 3D spatially flat universe in a 
totally relaxed state. 



TABLE I: Mass and energy fractions at a redshift of z = 1 
for flat and curved models and various values of fl™ 



Density 


flat 


curved 


flat 


curved 


flat 


curved 




0.02 


0.02 


0.03 


0.03 


0.04 


0.04 


Ha 


0.905 


0.899 


0.914 


0.905 


0.922 


0.912 




0.16 


0.104 


0.24 


0.162 


0.32 


0.225 




1.065 


1.003 


1.154 


1.067 


1.242 


1.137 



VII. CONCLUSION 

The 5D brane world of Moshe Carmeli has been has 
been applied to the accelerating expanding universe and 
the magnitude-redshift distance relation has been applied 
to the distance modulus data from the type la SNe mea- 
surements. It has been found that considering only the 
evolution of baryonic matter density as a function of red- 
shift, the resulting distance-redshift relation fits the data 
of the high-z supernova teams without the need for any 
dark matter. 

Astronomers and cosmologists have wondered whether 
the expansion of the universe is speed-limited in the usual 
Doppler sense. Here it has also been shown that the ex- 
pansion is, in fact, speed-limited. This means that even 
though the galaxies are assumed fixed within the expand- 
ing space, energy cannot be transported faster than the 
speed of light in vacuum, i.e. photons arriving at an 
Earth detector arc relativistically Doppler redshifted. 

Considering the range set on the present epoch bary- 
onic matter density of 0.007 < Q b < 0.041 8] with a 
best guess of Of, w 0.021 and the fact that necessarily 
fi m > 0, the best estimates of the matter density (which 
is only baryonic) from this analysis are O m = 0.0211q'q^4 
with the flat spacevelocity model. A similar result was 
obtained fitting a curved spacevelocity model. 

Even though the statistically derived standard errors 
are large, it is only necessary that f2 m be within the range 
of the locally measured baryonic matter density for the 
theory to fit the data. Therefore dark matter can play no 
part in Carmeli's description of the large scale structure 
of the universe. 

Even though it does not explicitly appear in Carmeli's 
spacevelocity metric, the vacuum energy contribution to 
gravity f2 A is a property of the metric and tends to unity 
as a function of decreasing redshift. This indicates that 
the universe, although always open because Q m < 1, is 
asymptotically expanding towards a spatially flat state, 
with a value of Qa + ~ 1-021 at the present epoch. 



APPENDIX 

This appendix deals with a more rigorous approach to 
the effect that curved spacevelocity has on matter density 
as a function of redshift. However this approach results 
in a transcendental equation, which could not be used 
in analytical curve-fitting. Nevertheless the analysis here 
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FIG. 7: (Color online) Mass density f2 as a function of red- 
shift z for the two models. Curve 2 (blue solid line) is Q for 
the flat spacevelocity model with Q m — 0.021. Curve 1 (red 
dashed line) is Q for the curved spacevelocity model derived 
in the Appendix with O m = 0.027 



validates the methods used in sections lllll I1VI 

If we start with the z-derivative of f2 from the right- 
hand side of lfTT|) . we get 



dQ = 30 m (l + zfdz. 
By the chain rule ljA.l|l can be written as 

dz 

dn = 3Q m (l + z) 2 — dr. 

dr 



(A.l) 



(A.2) 



follows from (|19|) that dz/dr = 1/ct. Now using l|19|) . in 
general, 

e £/2 

dfl = 3fi m (l + z) 2 dr (A.3) 

cr 

which when substituting e^/ 2 from (JSJ and (J7J becomes, 
(1 + z) 2 



= 3fL, 



2 Vcr 



(A.4) 



Finally by substituting r/cr = z, which is the Hubble 
law to the lowest order in z in the limit of no gravity, we 
arrive at an expression which is integrable and describes 
the matter density in the presence of curved spacevelocity. 



dfl ~ 3fL, 



y/l + (l-il)z 2 



-dz. 



(A.5) 



From (|A.5(I it is clear we recover JA.lt when il = 1 . So 
all that needs to be done is integrate l|A.5|) to get as 
a function of z. But it is not simply solved. It appears 
difficult to separate f2 in the integrand on the right-hand 
side from z. So assuming that in (|A.5(1 the value of 
is approximately constant or only slowly changing as a 
function of z we integrate and impose the boundary con- 
ditions that when f2 = 1 and z<l then = f2 m (l + z) 3 . 
To satisfy both conditions in IjA.lfl the limits of integra- 
tion must be from —1 to z. That is 



^2 1 - " 1 -' 3 t. 



(l + z>) 



>\2 



-i ^i + (i-n)z' 2 



-.dz'. 



(A.6) 



In a flat spacevelocity where Q — 1 and e^/ 2 = 1 it This results in 



£7 m 
2(1 -ft) 3 / 2 



3VT~£i ((4 + z)y/l + (1 - tyz 2 - 3V2" 



fij + 3(1 — 20) \ arcsinh\/l — fl + arcsinh(z\/l — Q) 

, (All 



from which it may be determined that O — ► f2 m (l + z) 3 
for arbitary z in the limit where il — > 1. 

In the limit where z — > 0, £1 — > f2 m (l + z) 3 where 
-> 1 and n -> ft m (0.958 + 3z) » O m (l + z) 3 where 
S7 — ?• 0. It is expected that the mass density can alway 
be approximated to yield Euclidean space locally. From 
the former, the boundary conditions are almost met and 
hence the approximation is reasonably valid. 

The solution to l|A.7|) is not analytical and must be 
solved numerically. The form of the transendental equa- 
tion does not lend itself to curve fitting as has been done 
in sections lllll IIVI However, using the software package 



Mathematica, it is possible to create a function that is the 
solution to such an equation as l|A.7|) . This was done and 
compared with the flat spacevelocity model. The result is 
almost identical to i|ll|) for redshift z < 2, that is, when 
Q m = 0.021 in the flat spacevelocity model, fi TO = 0.027 
in the resulting curved spacevelocity model. See figure 
where the density ft has been plotted for both models up 
to z = 2. 

This then means the matter density ft(z) in the curved 
spacevelocity model has approximately the same func- 
tional dependence on redshift z used in the main text. 
Therefore we are justified in using the flat spacevelocity 
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model or the approximate curved spacevelocity model. 
Only where z > 2 does the curve derived from (|A.7|I de- 
part from the flat space model of (|11[) . Also for z>2 the 
assumption, that be at most a slowly changing func- 
tion, breaks down and thus the integration in (|A.G!l is 



invalid. 

At these higher redshifts the unapproximated form of 
(I A. 4(1 must be solved where r/cr is substituted from lfTT)|) 
and v/c -> ? = ((1 + z) 2 - 1)/((1 + z) 2 + 1). But that is 
beyond the scope of this paper. 
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